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Abstract 

Dengue is one of the most important infectious diseases in the world, in terms of death and economic cost. Hence, the modeling 
of dengue is of great importance to help us understand the dynamics disease, and interfering with its spreading mathematical by 
the proposition of control methods. In this work, control strategies in an attempt to eliminate the Aedes aegypti mosquito, as 
well as proposals for the vaccination campaign are evaluated. In our mathematical model, the mechanical control is accomplished 
through the environmental support capacity affected by a discrete function that represents removal of breeding. Chemical control is 
carried out from the use of insecticide and larvicidal. The efficacious of vaccination is studied through the removal of a fraction of 
individuals, proportional to the rate of vaccination, from the susceptible compartment and its transfer to the recovered compartment. 
Our major find is that the dengue fever epidemic is only eradicated with the use of an immunizing vaccine because control measures 
directed against the vector are not enough to halt disease spreading. Even where the infected mosquitoes are eliminated from the 
system, susceptible mosquitoes are still present, and infected humans cause dengue fever to reappear in the human population. 
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1. Introduction 

In recent years the field of public health has been benefit from 
the use of mathematical models to the study of spread epidemic 
of infectious disease. One of the first mathematical model epi¬ 
demiology has been proposed by W. H. Hamer in 1906, who 
postulated that the development of an epidemic is related to the 
rate of contact between infected and susceptible individuals (H 
and currently this postulate is known as the principle of mass 
action O. 

In 1911 Ross has proposed one of the first models used 
for epidemiological study of malaria In 1957 Macdonald 
brought a greater biological realism to the Ross model, improv¬ 
ing the accuracy and interpretation of the model parameters. 
The contributions of Macdonald resulted in the creation of the 
model known as the Ross-Macdonald model n, in which two 
differential equations represent the dynamics of the disease re¬ 
gardless of the age structure of populations and the incubation 
period. 

In the cases of infectious diseases, situations may exist in 
which infected hosts once recovered again become susceptible 
to new infections or may become immune. However, there may 
be situations where the host contract the disease and do not ac¬ 
quire the disease elicited by the same agent again. In the first 
case, these diseases are usually caused by bacterial or protozoan 
agents Eia, such as malaria for example. The second case is 
due to viral diseases such as dengue Q. 

Currently, there is great interest in studying epidemic propa¬ 
gation, as in the case of dengue. This disease is considered the 
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most significant in terms of death and economic cost, since it is 
responsible for more deaths and infections in humans than any 
other arthropod-borne viruses ii. This scenario can be even 
worse, since the exact number of cases of the disease is not 
known. Several cases go unreported because they may clini¬ 
cally manifest with few or no symptoms a. Due to the large 
number of people exposed, and absence of medicines able to 
cure infected people, the increase of dengue epidemics have 
motivated an intensive study of mechanical and chemical con¬ 
trol mechanisms to used on eliminate the mosquito Aedes ae- 
gypti. 

The mechanical control, relying on educational campaigns 
for the removal of containers that accumulate water, results in 
a significant decrease in population size of the vector when 
this type of control is performed on any of the seasons ifTOl . 
Moreover, the mechanical control is effective only upon the 
awareness of the human population to identify and eliminate 
outbreaks of dengue in their own homes fTTl . Chemical con¬ 
trol, based on the use of insecticide and larvicide, was tested 
for the reduction of vector population |[II]|T2l. More recently, 
in parallel, vaccination strategies have been proposed as vac¬ 
cines. Studies have been done to develop a vaccine, evaluate 
and investigate its efficacy |[T3l[T^fT5]| . 

So, in this paper, we evaluate the set of these three strategies 
- mechanical control, chemical control, and vaccination - in or¬ 
der to assess the best way to eradicate dengue. Our study is 
based on the dengue epidemiological model proposed by Yang 
and Ferreira ca. We will model the effects by the use of me¬ 
chanical control mechanisms (parameters that decrease the rate 
of oviposition vector, known as carrying capacity) and chemical 
control (insecticide and larvicide). Besides the mortality rates 
of the vector we have also included in our study the spreading of 
the dengue virus in humans and suggest vaccination strategies 
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in a human population. 

The paper is organized as follows. Section Q introduces 
the ordinary differential equations of the model. In Section ([^, 
mechanical and chemical control strategies, as well as the fac¬ 
tors that justify the use of information based strategies in cli¬ 
mate and rainfall modelling are described. It is also highlighted 
the effects of vaccination campaign and reported a measure of 
the efficiency of these mechanisms. In section Q, the model 
solutions obtained by numerical integration are presented and 
the efficiency of each strategy is evaluate. Finally, in Section 
some conclusions are drawn. 


2. Mathematical model for dengue 

In the study of some epidemics, it is useful to separate the 
vector and host populations in some classes according to their 
stage in the infection cycle. In the case of dengue, for instance, 
it is useful to consider the existence of three classes for the vec¬ 
tor and four for the host. Such separation is suggested, because 
both the vector and the host contract the virus, go through a pe¬ 
riod of latency in which the individual has the virus but is not 
capable of transmitting it, became infective and, for humans, 
recovered from the disease uni. 

The dengue mosquito, the Aedes aegypti female, can be 
infected with the virus when feeding infected human blood. 
Once it have acquired the virus the mosquito, remains infected 
throughout its life. So it is feasible to divide the population of 
adult midges into susceptible, exposed and infected. Classes in 
turn, the human population can recover from the illness, and 
for that reason, in addition to the existing classes for the vector, 
there is the class of recovered humans d- 

The model in which the individuals are organized in these 
different classes from which they can entry or exit is known 
as 5'£’//?(susceptible - exposed - infectious - recovered) f\M . 
In this type of approach, a population is represented by state 
variables (subpopulation densities) associated to the different 
model compartments, which envolves on time according the 
rules of interaction. Differential equations can be used to model 
the rules, and give every subpopulation growth rates 12011 ^ . 

Our mathematical model, based on the model developed by 
Yang and Ferreira d, comprises 10 nonlinear first order ordi¬ 
nary differential equations. These equations include variables 
describing two classes of individuals: humans and mosquitoes. 
The phase of the mosquito in life cycle are represented in the 
system Q as follows: egg (E), larva (L), pupa (P) and adult 
(W), the latter being composed of three subpopulations differ¬ 
ing in its stage of infection, namely, susceptible (ITi), exposed 
(W 2 ) and infected (IT 3 ). 

The parameters of the model are related to biotic and abi¬ 
otic factors. They represent averages of values that takes into 
account all stages of vector in development, effective rates of 
natural mortality, and contact rates between the stages of infec¬ 
tion of the human population and the mosquitoes. 

The transition rates between the phases of egg, larva, pupa 
and adult are represented by the parameters ctj, (j = ej,p,w). 
In addition, the system 0 contains other rates: oviposition 


( 0 ), mortality of eggs, larvae, pupae and adult midges pt (i = 
/, p, w); which are functions of temperature, humidity, and rain 
precipitation. Thus these rates reflect abiotic factors. However, 
in the present analysis, these functions will be replaced by, con¬ 
stant average values, listed in table 0 . 
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Other parameters correspond to biotic factors, i.e. those that 
have direct interference of the human activity: carrying capac¬ 
ity, represented by the letter C' and the chemical control repre¬ 
sented by the function p'.. The carrying capacity can be mod¬ 
ulated through the elimination of breeding grounds as a result 
of educational campaigns to increase the population awareness, 
concerning the necessities to eliminate the aquatic population 
reservoirs. The chemical control is another parameter which 
takes into account the work of health agents responsible for 
preventing vector proliferation in every stage of its life cycle, 
from vector aquatic phase to adulthood. The population adult 
mosquitoes is composed only by females feeding blood. 

The equations for the humans population are coupled and 
separated in susceptible (s), exposed (e), infected (/) and re¬ 
covered (r) compartiments. In our model we take into consid¬ 
eration the presence of only one virus serotype. System 0 in¬ 
troduces a set of four first order, nonlinear ordinary differential 
equations for the human population densities. 
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3. Control strategies 

3.1. The influence of climate and temperature 

The climate is one of the factors that affects the environment 
where the Aedes aegypti is present Ea. The Aedes is found 
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primarily in tropical and subtropical areas in zones having rel¬ 
atively constant temperatures around 24^ C (231. The tempera¬ 
ture has direct and indirect influences on insect’s development 
and feeding (241 . Hence, the effects of temperature on the life 
cycle of insects read be modeled to understand the population 
dynamics of dengue vectors and can be exploited to develop 
appropriate vector control strategies. 

In Brazil, the fall season begins in March, winter in June, the 
spring starts in September, and the summer in December. How¬ 
ever, these four seasons are neatly marked only in the southern 
region, Sao Paulo, Mato Grosso do Sul, and in the mountain 
regions of Minas Gerais and Rio de Janeiro, occupying approx¬ 
imately 15% of the brazilian territory. In such regions sum¬ 
mer and winter temperatures are well-defined. In contrast, in 
the Amazonia, for example, does not exist significant variation 
in temperature and rainfall during the annual period, so basi¬ 
cally there are no seasons. In other regions, there are only two 
regimes: the rainy and the dry seasons. 

Taking into consideration the climatic variability in Brazil, 
we have divided a year into three periods, defined in terms of 
the maximum and minimum rainfall a major determinant of the 
mosquito population namely, favorable, unfavorable and inter¬ 
mediate. Indeed, the Aedes aegypti has a weak resistant at to the 
low temperature. Temperatures between 22^ C and 30^ C are 
favorable in its life cycle, and the extremes, below 18^ C and 
above 34^ C have negative effects on their development fISl . 

The most favorable period for the development of the vec¬ 
tor are the months of December, January, February and March, 
due to high temperatures and high rainfall. The colder months 
(June, July, August and September), are unfavourable for the 
development of Aedes aegypti. The intermediate period in¬ 
cludes the months of April, May, October and November (26l. 
The values that make up the template have daily values accord¬ 
ing to table Q and were taken from studies by Yang and Luz 

CSIIIII. 

3.2. Mechanical Control 

It is widely known that a population can not grow indefi¬ 
nitely because the environment has a limited capacity of re¬ 
sources to sustain their individuals. This ability is known as 
environmental supportability There population size of the 
dengue mosquito increases considerably rainy season due to the 
increase in the number of containers with standing water, ideal 
places for mosquito oviposition (28]| . Thus, the environmental 
support capacity of the dengue vector is primarily related to the 
amount of available sites for laying eggs, since the resources 
necessary for the maintenance of adults vectors such, as food 
and shelter, are virtually unlimited. 

In the model, described by equations Q, the environmen¬ 
tal support capacity is given by C' for the population of eggs, 
and is interpreted as the physical space available for mosquitoes 
oviposition. A method to control the growth of mosquito pop¬ 
ulation is to decrease C through the elimination of breeding 
grounds. This is usually done mechanically, at different times 
of the year, accordingly the climatic conditions. 

In order to take ino account this control strategy, our model 
assumes that the support capacity C is a function of a random 


variable C/. The effectiveness reate of the control performed at 
given moment. Specifically, 

c' = CiCfi^o (3) 

where Cis the supportability of the containers in the absence 
of control mesures and 

0 < G < 1 (4) 

varies according to the time of year. 

When Ci receives values close to 0, leading to a very small 
C', there are few locations for oviposition. In contrast, when 
Ci ~ 1, the support capacity remains virtually unchanged, and 
there is no reduction of mosquito breeding sites. In turn, the 
value of Cfixo varies with the annual seasons. 

Within the favorable period, characterized by high tempera¬ 
tures and humidity (29]| . a value Cpxo = 700 was fixed. Within 
(cold and arid) the unfavorable period fSUl we set Cpxo = 300. 
For the months that comprise the mild period, a supportability 
Cfixo = 500 is assumed. 

3.3. Chemical Control 

There are two types of chemical that hill the dengue vector: 
insecticides, that affect adult mosquitoes, and larvicides, that 
kill vector larvae. We assume that the chemical control does 
not act on the population of pupae (yu^ = 0). The insecticides is 
the most used instruments for vector control. In turn, larvicides 
are complementarly imployed to this control. Hence, in our 
model, these chemical methods are represented by the variable 
p'^ an p'l that affect winged mosquitoes and larvae, respectively. 
These mortality rates are exponential decay functions. They 
change with the period of the year as shown in table It is 
assumed that the chemicals are released during the first weeks 
of every month, and persists for 15 days in the environment. 
Also, at the first day the chemical release is the maximum, and 
at most 10% of the total amount of released chemicals remain 
in the environment. 

3.4. Vaccination Campaign 

Currently the fight against dengue is restricted to eradication 
methods of its vector. The development of new tools to com¬ 
bat dengue, as for example, the discovery of a vaccine requires 
more knowledge about the biological characteristics of the virus 
and its interaction with the hosts. 

Our model, described by equations [^and|^ does not take into 
account both people who could die by the disease and does not 
describe the existence of a vaccine that makes humans immune 
to dengue virus. 

The development of a tetravalent vaccine against the dengue 
virus is a promise of the World Health Organization (WHO), 
but efforts to develop vaccines have been hampered by the lack 
of an appropriate animal test subject (3T1l . 

In order to take into account these two features of the epi¬ 
demics, we extended our model as follows. The death of in¬ 
fected human introduced through the pd = 0.04 addition of a 
mortality rate in the infected compartment: 

d 

— i{t) = jh e{t) - (o-h + /i/z + Pd) Kt) (5) 

dt 
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Meaning of the Parameters 

Symbol 

Value (day ^) 

Oviposition 

0 

1 

Transformation from egg to larva(F) 

O-e 

0.330 

Transformation from egg to larva (U) 

CTe 

0.300 

Transformation from egg to larva (I) 

O-e 

0.200 

Transformation from larva to pupa (F) 

cri 

0.140 

Transformation from larva to pupa (U) 

o-i 

0.125 

Transformation from larva to pupa (I) 

o-i 

0.066 

Transformation from pupa to mosquito (F) 

o-p 

0.346 

Transformation from pupa to mosquito (U) 

CTp 

0.323 

Transformation from pupa to mosquito (I) 

O-p 

0.091 

Mortality rates of egg 

Ee 

0.050 

Mortality rates of larva 

El 

0.050 

Mortality rates of pupa 

Ep 

0.0167 

Mortality rates of mosquito (F) 

Ew 

0.042 

Mortality rates of mosquito (U) 

Ew 

0.040 

Mortality rates of mosquito (I) 

Ew 

0.059 

Force of infection human infectious ^mosquito susceptible 

Pw 

0.750 

Force of infection human susceptible^mosquito infectious 

Ph 

0.375 

Incubation period of virus (mosquito) 

Tw 

0.200 

Incubation period of virus (in humans) 

7h 

0.100 

Infectious (or recovery) period (in humans) 

o-p 

0.143 

Rates per-capita mortality human 

Eh 

0.000042 

Rates per-capita natality human 

En 

0.00042 


Table 1: The parameters and their values used in the model. The transformation rates in the aquatic phase as well as rate mosquitoes mortality are subdivided into 
periods: Favorable (F), Unfavorable (U) and intermediate (/). The values have units of day“^ 


Period 

Insecticide function 

Larvicide function 

Unfavorable 

Intermediate 

Favorable 

= 0.941 e-uouiw* 

= 0.960 
= 0.958 

= 0.884 

yu;(0 = 0.825 e-oooi4u 

yu;(0 = 0.810 


Table 2: Continuous functions representing the control by insecticide and larvicidal, at the favorable, unfavorable and intermediate periods. 


To represent a possible vaccination campaign, a fraction pro¬ 
portional to a vaccination rate v = 0.2 of susceptible individuals 
is removed from and added directly to the recovered compart¬ 
ment. Also, the efficiency of the vaccine is takes into account 
through a variable If ^ = 0, the vaccine is not effective and 
does not change the system, whereas, if ^ = 1, the vaccine is 
fully effective, and all the vaccined individuals became immu¬ 
nized. So, the effect of a vaccination campaing on the suscepti¬ 
ble and recovered populations are given by: 


— s(t) = jlh- 
dt 
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The duration of the campaign is varied to find out the time 
needed to eliminate the virus using a vaccination campaign. 


3.5. Efficiency 

For each control methods, its percentual efficiency E can be 
defined as 

l-^|100% (8) 

I 



where Ac represents the area under the affected mosquito curve 
population as a function of time when the controlmethod is used 
and As is the area under the affected total population curve 
without the use of the method. 


4. Results 

The results refer to the numerical simulation for a time pe¬ 
riod of approximately 20 to 50 years. The parameter values 
used by the table Q were adopted in order to describe coexis¬ 
tence between populations of humans and mosquitoes (aquatic 
and adult). Initially, we consider that the human population 
grows exponentially, as suggested by the United Nations (UN). 
Finally, simulation for a stationary human population is consid¬ 
ered. 

4.1. Human population with exponential growth 

Dengue is a disease recognized in more than 100 countries 
around the world Ea. According to the UN, over 70% of the 
world’s population live in underdeveloped countries. Asia, with 
a population of 4.1 billion inhabitants, and Africa, with 1.031 
billion inhabitants exhibit, both a high population growth rate. 
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which can be characterized as exponential 13^ . Since in those 
continents there is the presence of dengue viruses, there is much 
current interest in analyse the dengue epidemics in such expo¬ 
nentially growing populations. 


4.1.1. Mechanical control 

The Aedes aegypti female lays its eggs in places contain¬ 
ing stagnant water which nurtures the development of imma¬ 
ture stages of the vector. Consequently, if breeding grounds are 
eliminated, a decrease in the population size of the vector and, 
consequently, in cases of dengue fever will expected. In order 
to test the impact of the mechanical control on the population 
of adult vectors, a periodical decrease has been made in the car¬ 
rying capacity, following the rule described in section (3.2). 

Figure Q compares the results when the mechanical control 
was accomplished distinct over annual periods (favorable, un¬ 
favorable, intermediary and throughout the year) (deshed line) 
with those in which no control was performed (full line). The 
greatest efficiency was obtained when the mechanical control 
was carried out throughout the year, as shown ([2- J). Figure 
Q describes the evolution of the infected mosquito and human 
populations, when no control is performed (full line) and when 
the mechanical control is performed only in the intermediate 
period (dashed line). In the time interval between 240 and 360 
days, there was an increase in the number of disease case when 
compared to simulations in which there were no mechanical 
control. However, in the next interval (between 360 and 480 
days) this number was smaller. This variation in the number of 
cases is due to the interference of mechanical control. Reflexes 
are observed in the period after the intervention, since this type 
of control eliminates the aquatic stages of the mosquito. 


- Mechanical Control 

- Without Control 




Time in Days 




Time in Days 


Figure 1: Comparison between the population sizes of the dengue vector with 
and without mechanical control. The control is held only in (a) the most favor¬ 
able period, (b) unfavorable period, (c) intermediate period and (d) throughout 
the year. The case when mechanical control is performed (deshed line) is com¬ 
pared to the case in which no control type was performed (full line). 


When the control is performed only in one of the annual 
periods, an increase in the number of humans infected in pe¬ 
riod, even with a reduction in the number of infected vectors as 


shown in flgurej^ 

We can realize greater efficiencies through the use of me¬ 
chanical control throughout the year. However, the Aedes ae¬ 
gypti mosquito is not eradicated only using the mechanical con¬ 
trol and, consequently, the virus remains in the human popula¬ 
tion. Hence, the use of insecticide and larvicidal was tested in 
the following section in an attempt to exterminate the disease. 
The corresponding results are discussed below. 


- Without Control 

- Mechanical Control 




Figure 2: (a)Population size of mosquitoes infected with and without the me¬ 
chanical control in the intermediate period, (b) Population of humans infected 
when no control is performed and when the control is performed in the inter¬ 
mediate period. 


4.1.2. Chemical Control 

Christophers found that the mosquito in larval and adult 
stages are billed by of chemicals applied during their develop¬ 
ment El. We use the functions in the table (|^, in order to eval¬ 
uate the impact of insecticides and larvicidaes on the control of 
dengue epidemic. In addition, to discriminate the impact of dif¬ 
ferent types of chemical control, simulations were performed 
considering only one type of control at a time. 

Figure 0 presents the population size of the Aedes ae¬ 
gypti mosquito when an insecticide was applied. As can be 
seem, the use of the insecticide was not effective in eliminating 
the mosquito population, and when the application was made 
throughout the year (Figure 3-d), the impact on the total popula¬ 
tion was greater. However this control had only 21% efficiency, 
as calculated according to the equation 

In (|^ are shown the results for the application of larvicidal. 
The approximate efficiency is 5%, and the population practi¬ 
cally evolves as if no control mechanism was used. This indi¬ 
cates that such control is not feasible in an attempt to eradicate 
the Aedes aegypti mosquito. 

Figure 0 shows the effect of chemical control on the pop¬ 
ulation of infected humans. The result is independent on the 
chemical applications. Even with a decrease in the number of 
individuals at a given time, this number increases again latter 
leading do a control efficiency of less than 1%. This may indi¬ 
cate that the density of infected mosquitoes does not interfere 
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Figure 3: Mosquito population sizes after applying the insecticide only in (a) 
the most favorable, (b) unfavorable, (c) intermediate period, and (d) throughout 
the year. 
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Figure 5: Comparison of the population size of humans infected when no con¬ 
trol is carried out and when the chemical and mechanical control are performed 
together throughout the year. 
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Figure 4: Comparison of the population size of mosquitoes after application of 
larvicidal when the control is held (a) only in the most favorable period, (b) the 
unfavorable period, (c) the intermediate period and (d) throughout the year. 


with the virus permanence into the system, but simply there 
is a critical population size that guarantees virus support. Ac¬ 
cordingly, in countries in an endemic state, most of their dis¬ 
tricts presented a vector index below the value required for the 
risk of dengue fever epidemics 1291 . The incidence of infection 
was high even in areas with low numbers of positive dengue 
mosquito breeding sites ESI. Furthermore, small breeding are 
as with mosquito outbreaks in their neighborhoods can ensure 
the presence of dengue virus in all other districts of a city Oil. 

The long-term control done by chemical, such as those con¬ 
sidered here, should be analysed with caution because its con¬ 
tinued use can generate resistant vectors populations urn . In 
addition, according to the technical note n" 109/2010 1 ^ . it is 
important to consider that there are no insecticides suitable for 
use in public health, since all these products are used specifi- 


4.1.3. Mechanical and chemical control 

Since the use of insecticide and larvicidal separately was in¬ 
efficient to eliminate the population of dengue adult vectors, 
we considered a scenario in which the two chemicals were ap¬ 
plied simultaneously with the mechanical control throughout 
the year. 

The results are shown in fig. 0 for different control strate¬ 
gies. When chemical control was combined with the mechani¬ 
cal control, the efficiency is 86% (blue curve). In comparison, 
if only the mechanical control (green curve) was used the effi¬ 
ciency fell to 26%. However, even if both types of control are 
used, it has not been possible to eliminate dengue from the sys¬ 
tem, still leaving the human population at risk of contracting 
the disease. 



Figure 6: (Color online) Comparison between the different applications of 
chemical and mechanical controls. The combined use of chemical mechani¬ 
cal controls resulted in an 86% efficiency (blue curve). In turn, the efficiency 
fell to 26% it only the mechanical control (green curve) is used. 


4.1.4. Death by disease and vaccination campaign 

In previous simulations the fact that people can die due dis¬ 
ease complications was not considered. Dengue fever, espe¬ 
cially in its most severe form, is hemorrhagic, and may have a 
lethality greater than 20% l(3^ . Hence, in order to increase the 
model’s realism, we add to the human infected compartment. 
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eq. 0, the rate of death by illness Since, the Epidemio¬ 
logical Bulletin of the Ministry of health in the year of 2005 
reported a 4% lethality level for dengue, we adapt jid = 0.04 [? 

]. 

This new factor has a quantitative impact on disease out¬ 
breaks, leading to outbreaks delays at long-term, as can be seen 
in figure 0. 
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Figure 7: Evolution in time of the human infected population with and without 
the disease mortality rate parameter jid- 

Because, the chemical an mechanical controls were unable to 
eliminate the dengue virus circulation in the human population, 
we tested the effects of an eventual vaccine against dengue fever 
on its spreading dynamics. In order, to simulate a vaccination, 
a fraction, proportional to the vaccination rate, v = 0.2, is with¬ 
drawn from the compartment of susceptibles, equation ([6k and 
directly added into the recovered compartment, equation^. 

It is widely known that vaccines licensed by the health agen¬ 
cies are not 100% effective. Due to this fact, the efficiency of 
the vaccine was also considered and represented by a variable 

^ = 0 corresponds to an ineffective vaccine that does not de¬ 
velop immunity, whereas ^ = 1 means that the vaccine is fully 
effective, immunizing all the individuals received the vaccine. 

With a vaccination rate fixed at 0.2, and an efficiency main¬ 
tained at ^ = 1 the time during which the campaign should 
be maintained in order to eliminate the virus was estimated. 
Figure shows that a campaign active for three consecutive 
months eradicated the virus from the population. Furthermore, 
figure (j^ shows the evolution in time of the population infected 
humans for different values of ^ for a campaign lasting three 
months. A vaccine that is 60% efficient does not eliminate the 
virus and with 70% there is a significant decrease in the number 
of cases over time. In order to eliminate the dengue fever of the 
human population, at least 80% efficiency is required. Larger ^ 
values, such as 90%, only differ in the time spent to eradicate 
the virus. 

4.2. Human population without exponential growth 

As variously mentioned, dengue is present in many countries, 
including Brazil and other South American countries, where 
the rate of population growth is not enough to reach levels as 



Figure 8: Evolution in time of dengue epidemic for different durations of the 
vaccination campaign. 
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Figure 9: Comparison of different rates of efficiency for the vaccine against the 
dengue virus: (a) describes the efficiency of 90%, (b) 80% efficiency, (c) 70% 
efficiency, and (d) 60% efficiency. 


high as that of the African countries. According to a study by 
the Brazilian Institute of geography and statistics (IBG^, the 
Brazilian population is decreasing its growth rate since the 60'5', 
when this rate was 3.04%. The IBGE estimates that this reduc¬ 
tion in the growth rate still should last for many years reaching 
the so-called “zero growth” in 2039 ll40l . Hence, we also ad¬ 
justed the birth rate of the human population to a value closer to 
the mortality rate. So, in these simulation the birth rate is now 
fixed in 0.000045. 

Considering this new birth rate for the human population and 
the use of mechanical control throughout the year, an efficiency 
of 80% was reached. However, as shown in figure ( [T0| ), the 
reduction of mosquitoes population does not decrease the in¬ 
fected human population. The behavior remains the same when 
the control was not used. We note that even with a low num¬ 
ber of infected mosquitoes, the dengue virus still remains in the 
human population. This is possible because, as proven through 
histological studies, the dengue mosquito feeds on blood sev¬ 
eral times during a single cycle consequently, a single female 


^IBGE: Brazilian population ageing in rhythm accelerated 
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can contaminate several humans in the same cycle 1 411 . 
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Figure 10: Comparison of mosquito (total and infected) and infected human 
infected population when no control is in effect and when it is done using the 
mechanical control throughout the year with birth rate jid = 0.000045. (a) 
The total mosquitoes population, (b) the population of infected mosquitoes, (c) 
population of infected humans. 


Since only the mechanical control was enable to eliminate 
the dengue fever epidemic, a chemical control was applied in 
order to check its efficiency taking in to account this new birth 
rate. The system response to the use of insecticide and larvi- 
cidal applied in different periods of the year is shown in figure 
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Figure 11: Evolution in time of mosquito and humans population with the use 
of chemical control and birth rate jid = 0.000045 for the favorable, unfavorable 
and intermediate periods: (a) is the total population of mosquitoes, (b) is the 
population of infected mosquitoes, (c) is population of infected humans. 

Now we can notice that the chemical control is not efficient 
in eradicate the dengue viruses regardless of the period in which 
it is being applied, eliminating all infected mosquito in the sys¬ 
tem, but this does not reduce the levels significantly below the 
susceptible vector’s population, as Luz ITTIl . In the work done 
by Burattini and et. al, it was possible to evaluate the use of 


chemical control, which had a considerable impact to prevent 
the spread of dengue fever epidemic ifT^ . 

But the persistence of Aedes aegypti in the system, even if 
it is not infected, the dengue fever may reappear at any time if 
any infected person is introduced, and so the dengue epidemic 
can re-emerge in the human population, since the vector of the 
virus is still present. 

It is known that the vector elimination campaigns by means 
of chemical agents often occurs in the favorable period of the 
year. Accordingly, we permed the chemical control in this pe¬ 
riod to see how the system respond to the emergence of infected 
people. 

Figure ( p^ shows the result of simulations in which 10 in¬ 
dividuals infected with the dengue virus appear at two different 
periods of the year. If these infected individuals come when 
chemical control is applied in the favorable period (continuous 
line), dengue does not reappear and the system remains free of 
the disease. On the contrary, if they are introduced in the un¬ 
favorable period, when there was no chemical control (dashed 
line), we have periodic outbreaks of the disease are observed. 
This suggests that a control campaign performed only in a pe¬ 
riod of the year may not be effective in eliminating dengue to¬ 
tally, since the virus can reappear through the emergence of in¬ 
fected people. 


- Favorable 

-Unfavorable 



Figure 12: Evolution in time of populations of infected (a) mosquitoes (b) 
humans for a scenary in with 10 infected humans appear in a system free of 
viruses. The analyses were made in the favorable and unfavorable periods. 

With the rate pd = 0.000045 the human population does 
not grow exponentially, probably because the elimination of 
infected persons through the use of chemical control becomes 
possible. But even with the complete elimination of the popula¬ 
tion of infected vectors, it has not been possible to eliminate all 
the susceptible mosquito population in the system as shown in 
figure The dengue fever may reappear at any time, unless 
there is a vaccine to protect the population against the virus. 

4.2.1. Vaccination campaign 

In this simulation the population does not exhibit a high 
growth rate and the vaccination campaign will be represented 
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by the withdrawal of a fixed number of susceptible individuals 
(ij/ = 200) by adding them into the recovered population. Thus, 
the vaccination rate is no longer proportional to the number of 
susceptible individuals, and is described by equations and 
( p^ . In this case, a deemed efficiency ^ = 80% for the vaccine 
is considered since it is the immunization threshold acceptable 
by health agencies. 


d s{t) 
dt 




(9) 


d r{t) 
dt 




( 10 ) 


Figure ( p^ shows how the system behaves under immuniza¬ 
tion of 200 people per day for 1 month, with different times 
of duration of the campaign. It can be observed that, if the 
campaign remains for only 5 consecutive years, dengue virus 
re-emerge in populations of mosquitoes and humans and the 
epidemic reappears. However, if the duration of the campaign 
is 10 years, it is possible to eradicate the dengue virus from 
the human population for at least 30 years, in a total simula¬ 
tions time fixed at 50 years. To find out the time required for 
the permanent eradication of the virus from the system, more 
simulations considering longer simulation times, as well as the 
possibility of the emergence of other virus serotypes must be 
carried out. 



Figure 13: The vaccination campaign to immunize 200 individuals per day with 
different periods of duration. 


Several studies show that a vaccine that protects the popula¬ 
tion from the dengue virus is about to be created. An attenuated 
tetravalent vaccine virus has been tested in order to character¬ 
ize the interactions between the different serotypes of dengue 
virus, concluding that it is possible to create an effective vaccine 
against all four serotypes of the dengue virus 1421 . However, 
an important point to be considered is whether live, attenuated 
virus vaccines do not offer risks to human health, since it is not 
clear that the presence of four types of strains into a single vac¬ 
cine might increase the risk of developing severe forms of the 
disease ia. 


5. Conclusions 

Mathematical models are quite useful for observing the dy¬ 
namics of infectious diseases such as the dengue fever. The cre¬ 
ation of mathematical models that simulate the interaction be¬ 
tween host-vector generates momentum in the system. Through 
the compartmental model proposed by Yang and Ferreira was 
used, it was possible to introduce small changes in order to test 
new control strategies. Testing the mechanical and chemical 
control verified the vaccination strategies for the human popu¬ 
lation. 

The mechanical and chemical control were effective in re¬ 
ducing the population of Aedes aegypti but assured that the dis¬ 
ease has been eliminated from the system. The control that used 
only the larvicidal had a very low contribution in this reduction, 
probably because the time interval between the implementation 
and its impact on the population of adults is very high, making 
this control to be less appropriate in controlling an epidemic of 
dengue fever. The mechanical control is much more efficient, 
which indicates that awareness campaigns are more useful than 
the use of controls through insecticide and larvicidal. 

When considering a human population birth rate equal to 
0.00042/it was possible to simulate an exponential growth 
dynamics of this population. However, this consideration does 
not affect the effectiveness of mechanical control. In this case 
the application of mechanical control must be constant through¬ 
out the year. Once this control is active only in time after his 
application, it generates an increase in the number of cases in 
the other two periods, as there is no control being carried out. 

With respect to the use of chemical control, there was an in¬ 
terference mosquito population in a yield of 21% for the insec¬ 
ticide and larvicide only 5%. When the chemical and mechani¬ 
cal control were used together, there was an efficiency of 26%. 
In the latter case, even with greater efficiency, the use of these 
controls did not influence the dynamics of dengue in the human 
population, who remained infected. Importantly, the chemicals 
used in the simulations is fictitious and quite specific. Sim¬ 
ulations with insecticide and larvicide with different duration 
should be conducted so that we can guarantee these results. 

When the use of both controls was compared to situation 
where no control was accomplished, it resulted in an efficiency 
of 86%, compared to the situation where only the mechanical 
control was carried out in which the efficiency value fell to 26%, 
thus showing that the mechanical control is more efficient in the 
elimination of mosquito population than the chemical control 
even when applied throughout the year. 

When the birth rate was changed to 0.000045/Jay, it was no 
longer possible to observe an exponential growth for the human 
population. In this case the efficiency of mechanical control 
remained at 80%. In the same way for the previous simula¬ 
tions, there was no significant reflection on human population 
infected. The chemical control can eliminate the infected pop¬ 
ulation but the outbreaks do not end with a susceptible popula¬ 
tion. With 10 individuals in the system infected with the dengue 
virus in the period in which the control is not being applied in 
the unfavorable time, periodic outbreaks of the disease occur. 

The only way to contain the dengue epidemic, regardless of 
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the form of population growth, is to create a vaccine to im¬ 
munize the population against the virus. The population grew 
exponentially during the immunization vaccination campaign 
for a fraction proportional to a vaccination rate of individuals, 
y = 0.2. For this campaign to become really effective it is nec¬ 
essary to have a duration of three months and an efficiency rate 
of at least 80%. 

It was possible to note the existence of annual cycles and a 
strong seasonality of epidemic cycles 1221 . A major concern for 
the safety of live attenuated virus vaccines is the genetic stabil¬ 
ity. The use of molecular tools is required to ensure quality, 
safety and authenticity of genetic vaccine (441 . 

Already when the exponential growth is removed, the cam¬ 
paign is based on daily immunization ij/ = 200, of people, with 
an efficiency of 80%. In this case it is necessary to maintain the 
campaign for at least 10 years to ensure at least 30 years of a 
system free of the disease. 

Variable jid was introduced representing the portion of in¬ 
fected who die from complications related to dengue fever. The 
addition of this new factor caused the outbreaks of the disease 
might arise later, but as the simulation time is fixed at 50 years 
this variable can have a major stake, because it may delay the 
outbreak so that it will not reappear during the time examined. 
However, field surveys are still necessary in order to better un¬ 
derstand the mechanisms responsible for generating these dy¬ 
namics. In addition, more simulations must be carried out in 
order to consider a longer time of observation, while factoring 
in the possibility of the emergence of other serotypes of the dis¬ 
ease and biological ways of eliminating the vector. 

Further studies are needed to understand the mechanisms re¬ 
sponsible for generating the dynamics of dengue in order to 
learn how to fight it. The need for a vaccine is essential for the 
elimination of dengue viruses in human populations, thus it is 
expected that studies about immunization are satisfactory and 
achieve a positive result as soon as possible. 
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